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ABSTRACT 

Ultraviolet imaging with the Galaxy Evolution Explorer (GALEX) has re- 
vealed an extensive sample of UV-bright stellar complexes in the extreme outer 
disk of M83, extending to about four times the radius where the majority of 
H II regions are detected (Rhii = 5.1' or 6.6 kpc). These sources are typically 
associated with large-scale filamentary H I structures in the warped outer disk 
of M83, and are distributed beyond the galactocentric radii at which molecular 
ISM has yet been detected. We present measured properties of these stellar com- 
plexes, including FUV and NUV magnitudes and local gas surface density. Only 
a subset of the outer disk UV sources have corresponding H II regions detected 
in Ha imaging, consistent with a sample of mixed age in which some sources are 
a few Myr old and others are much more evolved (~ 10 8 yr). 

Subject headings: galaxies: individual (M83) — galaxies: evolution — ultraviolet: 
galaxies 
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1. Introduction 

The Galaxy Evolution Explorer (GALEX), with its 1.25° field-of-view and sensitivity to 
stellar populations younger than a few hundred Myr, is remarkably well-suited for addressing 
the topic of star formation at large galactocentric radii. Recent star formation within such 
environments has been detected in deep Ha and broadband observations of a few galaxies 
(NGC 628, NGC1058, NGC6946: Ferguson et al. 1998a, also Lelievre & Roy 2000 for 
NGC628; M31: Cuillandre et al. 2001; NGC6822: de Blok & Walter 2003). One of the goals 
of the GALEX Nearby Galaxy Survey (NGS: Bianchi et al. 2004) is to conduct a census of 
outer disk star formation in a statistically significant sample of galaxies. 

The confirmed occurrence of star formation in the outer disk of ordinary spiral galaxies 
has far-reaching implications. First, if massive stars are forming at large radii, then the 
outer disk continues to undergo chemical enrichment. This may clarify the origin of dust 
detected in the extended HI disks of spiral galaxies (Popescu & Tuffs 2003, A&A 410L, 
21P). Secondly, radiative and mechanical feedback processes associated with massive stars 
are expected to be more efficient in the comparatively rarified interstellar medium (ISM) and 
uncrowded environment at large radii. In particular, the unresolved issue of H I production 
via H 2 dissociation (Smith et al. 2000, Allen 1986) can be studied in detail. 

The presence of recently formed stellar complexes at large galactocentric radii also pro- 
vides a simplified laboratory for investigating the star formation threshold, the (minimum) 
gas column density required for star formation to occur spontaneously. The precise workings 
of this mechanism remain a matter of debate. Kennicutt (1989) demonstrated that star for- 
mation is rare in disk environments having total gas surface density £ 9as < 5 — 10 M Q pc -2 . 
Martin & Kennicutt (2001) showed the observed threshold gas density varies by about an 
order of magnitude at the "edge" of the star-forming disk (R = Rhii) in different galaxies, 
but the ratio of to a critical density (which depends on the physical and dynamical 
conditions of each galaxy in particular) is approximately constant at this position. 

Deep Ha imaging by Ferguson et al. (1998a) and Lelievre & Roy (2000) showed that 
star- formation regions only a few Myr old exist beyond two optical radii (R25) in NGG 628, 
NGC 1058, and NGC 6946. Lelievre & Roy (2000) demonstrated that the H II region 
luminosity function (LF) slope is significantly steeper in the outer disk, relative to a galaxy- 
wide sample, possibly indicating a genuine change in the initial cluster mass function (CMF) 
or a population which is preferentially older. Ferguson et al. (1998a) noted that the outer 
H II regions appear small, isolated, and faint compared to inner disk analogs. 

In this Letter, we present far-UV (FUV) and near-UV (NUV) wide-field GALEX im- 
agery which reveals an extensive population of recent star-forming regions in the extreme 
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outer disk of M83, out to ~ 3x the Holmberg radius. Some of these sources were discovered 
by Wamsteker et al. (1983) in deep red (127-04+RG630) plates, but observed properties 
were never published. We analyze the GALEX data in conjuction with visible and A21cm 
imagery, to place the outer disk sources in context relative to the underlying stellar pop- 
ulation and extended H I disk. The GALEX data is of particular importance because it 
demonstrates that Ha observations will fail to detect a significant population of moderate age 
clusters in the outermost disk of spiral galaxies, leading to an underestimate of galactic evo- 
lution in these locales. Processes such as chemical enrichment (Ferguson et al. 1998b), disk 
heating (Sellwood & Balbus 1999), and ISM phase balance (Allen 1986) may be significantly 
influenced, given that the (predominantly) B star population traced by GALEX effectively 
drives these feedback mechanisms. GALEX NGS observations show that M83 is not unique 
in terms of outer disk UV morphology. We will analyze UV-visible-IR outer disk properties 
of a representative galaxy sample in a subsequent paper. 



2. Observations and Data Analysis 

M83 (NGC5236) was observed with GALEX on 7 June 2003 for a single orbital night as 
part of the NGS, with an exposure of 1352 seconds in both FUV (1350-1750 A) and NUV 
(1750-2750 A). Figure 1 shows the color-composite GALEX image of M83. Two bright stars 
eastward of M83 prevented us from centering the target in the field, but our data includes 
most of the outer disk. RMS sensitivity is 26.6(26.8) AB mag for FUV(NUV), or 27.5(27.6) 
AB mag arcsec -2 expressed in units of surface brightness evaluated at the scale of the PSF. 

Boissier et al. (2004, this volume) also analyze the GALEX M83 dataset, emphasizing 
the inner disk in regard to extinction and chemical evolution models. 

Courtesy of P. Crosthwaite, we gained access to the A21cm VLA observations discussed 
by Tilanus & Allen (1993). For the inner 14' x 14' of M83, we also used the CO map of 
Crosthwaite et al. (2002) to estimate S gas as a function of position. Figures 1, 4, and 5 
pertain to the H I and CO observations. Narrowband Ha and R-band imaging of M83 was 
obtained using the Swope 1-m telescope at Las Campanas Observatory (see Fig. 5) 

We adopt a distance of 4.5 Mpc to M83 (Thim et al. 2003), thus the ~ 5" GALEX PSF 
corresponds to a linear scale of 110 pc. Accordingly, GALEX detections in M83 are generally 
not single clusters, but rather aggregates possibly consisting of a few discrete (unresolved) 
clusters. 

Photometry of the discrete UV sources within M83 was accomplished using Source 
Extractor (Bertin & Arnouts 1996) to define an isophotal aperture for every object detected 
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within our NUV image. These irregular region boundaries were then applied to both the NUV 
and FUV data for the measurement of background-subtracted source flux, after applying a 
distortion correction to the FUV image to ensure registration with our NUV data. The 
background for every source was estimated locally (within ~ 250 pc). We measured median 
surface brightness as a function of galactocentric radius as described by Bianchi et al. (2004, 
this volume) with % = 25° and PA = 226° (Crosthwaite et al. 2002). 

3. Recent star formation in the 
outer disk of M83 

M83 is considered prototypical of those galaxies showing a well-defined edge of the star- 
forming disk (eg. Martin & Kennicutt 2001). The surface density of H II regions declines 
abruptly near Rhii = 5.1' (6.6 kpc) radius (cf. Ruoim = 7.3', 9.6 kpc), suggesting that 
massive star formation is now confined to M83's inner disk, despite the prominence of a 
rather extended H I disk (Rogstad et al. 1974, Huchtmeier & Bohnenstengel 1981, see also 
Fig. 1). Our GALEX observations modify this view. While the distribution of UV-bright 
stellar complexes is highly concentrated in the area occupied by H II regions (tracing O 
stars), many UV sources (tracing O and B stars) are easily recognizable far beyond the 
"edge" of this distribution. Furthermore, diffuse UV emission is also detected beyond the 
optically-bright star-forming disk. 

3.1. Outer disk, UV-bright stellar complexes 

Our GALEX observations of M83 revealed > 100 UV-bright star-forming regions and 
stellar clusters beyond Rhii- Figure 1 shows that the outer disk star-forming regions appear 
preferentially located on local maxima or filaments in the structure of the warped H I disk. 
The distribution of GALEX sources extends to the limit of the field surveyed at the VLA 
by Tilanus & Allen (1993). We suspect that additional M83 star-forming regions may be 
located in the remainder of the H I envelope (Huchtmeier & Bohnenstengel 1981), but this 
is challenging to confirm without corroborating H I morphology. 

Upon inspection of POSS2 images, we found that the majority of the outer disk GALEX 
detections could be discerned in IIIaJ+GG385 (blue) plates. However, the apparent contrast 
of the sources is lower in the visible than at UV wavelengths. An undetermined fraction of 
the outer disk star-forming complexes were discovered by Wamsteker et al. (1983) in their 
analysis of deep red (127-04+RG630) Schmidt plates. To our knowledge, this is the only 
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mention of the outer disk sources in this galaxy prior to the GALEX data. 

Figure 2 presents the FUV luminosity function for M83 star-forming complexes in var- 
ious disk locations. The outer disk complexes (dashed line) appear systematically fainter 
than their inner disk counterparts (dotted line). The median FUV luminosity of the inner 
disk population is a factor of 3 greater than the corresponding figure beyond Rhii (only 
correcting for foreground extinction). Internal extinction is likely more severe in the inner 
disk. Thus, the intrinsic luminosity difference may be even higher than that shown in Fig 
2, although we caution that blending and photometric incompleteness (more significant at 
R < Rhii) would have an opposite effect. The shape of the apparent FUV LF differs be- 
tween the inner and outer disk populations. Our M83 FUV luminosity function of outer disk 
stellar complexes is steeper than for the inner disk population. A similar trend was noted 
in NGC 628 for the Ha LF of H II regions (Levievre & Roy 2000). Because our GALEX 
measurements represent the aggregate flux of more than one independent cluster (or simply 
an extended distribution of massive stars), it is premature to interpret the precise slope 
of these FUV LFs in terms of the CMF. We note that blending tends to flatten observed 
LFs, but the magnitude of this effect may vary with environment. High-resolution follow-up 
observations would resolve ambiguities related to the FUV LF and median luminosity esti- 
mates. Nevertheless, the maximum UV luminosity exhibited by the population of outer disk 
sources implies an upper limit of about 10 5 M© for the most recently formed complexes (see 
also Fig. 3). 

We compare the observed (NUV, FUV-NUV) color magnitude diagram (CMD) for UV- 
bright complexes with population synthesis models of Bruzual & Chariot (2003) in Figure 3. 
M83's outer disk sources have FUV-NUV colors generally consistent with ages from a few 
Myr up to 400 Myr and are generally less massive than sources with R < Rhii- 

Only a subset of the outer disk UV-bright sources have detected counterparts in Ha 
imagery. This effect could be partially due to a combination of limited sensitivity and the 
lower expected emission measure of H II nebulae in a low-density environment, but supports 
the view that our GALEX observations trace a diverse population, comprised of both young 
and moderate-age complexes. A distinct advantage of UV observations is to directly trace 
stellar, rather than nebular, emission. 

Figure 4 presents the distribution of total gas surface density at the locations of all 
inner(outer) disk sources using a dotted(dashed) line. The local surface densities at which 
outer disk star-forming complexes appear to be forming are typically subcritical, with respect 
to the disk instability criterion (Martin & Kennicutt 2001). However, we caution that the 
outer disk measurements of E 9as are lower limits in two important ways. Most obviously, the 
estimate of surface density at positions beyond ~ Rhii is based only on H I observations, 
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neglecting the contribution from molecular gas traced by CO (unmapped outside of 14' x 14'). 
Also, the resolution of the H I aperture synthesis map is coarse in comparison to the scale of 
our detections (50" vs. 5"). Higher- resolution 21cm observations could plausibly boost the 
apparent gas surface density in the local environment of the UV-bright complexes. 

3.2. Underlying low surface brightness disk 

The discrete outer disk stellar complexes lie projected against a fainter, slowly-varying 
field population. To place the UV-bright complexes in the context of M83's overall star 
formation history, we have plotted the median FUV, NUV and Ha surface brightness (Fig. 
5), plus FUV-NUV (Fig. 6), as a function of galactocentric radius. The underlying stellar 
population becomes substantially bluer (in FUV-NUV) with increasing radius. Dust is 
known to exist in the outer disks (and even in the extended HI disks) of spiral galaxies 
(Popescu & Tuffs 2003) and thus is likely to be present in the outer disk of M83, too, 
particularly given the star formation traced there by GALEX. Accordingly, scattered light 
may contribute to the emission at each position (see also Popescu et al., this volume), 
tending to smooth the UV profile. Nevertheless, the FUV-NUV profile (Fig. 6) shows a 
sequence of breaks within the bright disk of M83. Outside of Rhii the FUV-NUV color 
becomes bluer almost monotonically, indicating recent star formation dominates the UV 
emission. The median Ha surface brightness profile (Fig 5.) is observed to drop off much 
more rapidly than either UV band. Given that the dynamical timescale beyond Rhii is > 60 
Myr, the dropout of Ha emission compared to UV cannot likely be attributed to a (recent) 
coherent decrease in the SFR at large radii. FUV-NUV measurements for discrete sources 
are overplotted on Fig. 6. Within 10' galactocentric radius, the UV-bright star-forming 
regions generally appear bluer (younger) than or consistent with the FUV-NUV color of the 
underlying diffuse emission. 

GALEX (Galaxy Evolution Explorer) is a NASA Small Explorer, launched in April 
2003. We gratefully acknowledge NASA's support for construction, operation, and science 
analysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes 
Spatiales of France and the Korean Ministry of Science and Technology. 
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Fig. 1. — GALEX FUV and NUV color-composite image of M83, highlighting the newly 
re-discovered sites of recent star-formation in the extreme outer disk. We show FUV in blue 
and NUV in red (along with their average in green). A red contour represents the extent 
of the H I distribution detected by Tilanus & Allen (1993). The blue H I contour is drawn 
at 1.8 x 10 20 cm -2 , corresponding to ~ O.4M pc -2 . We also include a yellow contour from 
Crosthwaite et al. (2002), indicating total neutral gas surface density, T, gas of 1OM pc~ 2 . 
The yellow S sas contour lies immediately within Rhii- The projected spatial extent of the 
field (44') is 57.6 kpc (Rhii = 6.6 kpc) at a distance of 4.5 Mpc. 
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Fig. 2. — GALEX FUV luminosity functions of the UV-bright stellar complexes in M83, 
having photometric error < 0.25 mag. We show the FUV LF for all sources with a solid 
line, inner disk (R < Rhii) sources with a dotted red line, and outer disk (R > Rhii) 
sources with a dashed blue line. Outer disk complexes appear typically fainter and have a 
steeper slope in the FUV LF, relative to inner disk objects, although detailed assessment of 
completeness and blending is needed to bolster these statements. 
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Fig. 3.— GALEX (FUV-NUV, NUV) color-magnitude diagram for M83's UV-bright stellar 
complexes (same sample and radial classification as Fig. 2). Inner disk sources (R < Rhii) 
are plotted with red circles, whereas outer disk sources (R > Rhii) are indicated with blue 
diamonds. Our data points have been corrected only for Galactic foreground extinction, 
E(B-V)=0.066. The black curves show model predictions for starburst populations of 10 3 , 
10 4 , 10 5 , and 10 6 M as a function of age (log age = 5.7, 6.0, 6.3, 6.5, 6.6, 6.75, 7.0, 7.3, 
7.6, 7.75, 7.9, 8.0, 8.7, 9.0 indicated for 10 6 M ). Dashed lines indicate the same models 
reddened with E(B-V)=0.5, a value slightly higher than the median E(B-V) determined for 
a large sample of clusters in M83 (Harris et al. 2001). 
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Fig. 4. — Local S gas at sites of recent star formation in M83, evaluated at 1.1 kpc (50") 
resolution. The overall distribution (solid) appears largely bimodal, with inner disk UV 
sources (within Rhii, dotted red) forming in locales with total surface density greater than 
~ 1OM pc~ 2 . Recently formed, outer disk stellar complexes (dashed blue) exist within 
environments characterized by local surface density about an order of magnitude smaller, 
although the CO map used to compute the H 2 contribution does not extend to large galac- 
tocentric radii. High-resolution H I and CO observations toward the outer UV-bright com- 
plexes are critically needed. Such data would enable investigation of H I production via H 2 
photodissociation (Smith et al. 2000), as well as a fundamental check on the star formation 
threshold (Kennicutt 1989). We reiterate that this plot shows local measurements of S 5as , 
whereas Fig. 5 presents azimuthally averaged surface density as a function of galactocentric 
radius. 
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Fig. 5. — Radial profiles of the median FUV (blue), NUV (red), and Ha (black) surface 
brightness in M83, plus average surface densities derived from H I (green) and CO (orange) 
data. Azimuthally averaged T> gas is shown with a dashed line. The UV profiles are presented 
in absolute units, while the Ha profile is arbitrarily normalized. We indicate Rhii with 
a vertical line, and the UV sky background with horizontal lines. Slight deviations from 
flatness across the GALEX field, coupled with the fact that M83 is positioned on the edge 
of the field, make it impractical to follow the profiles beyond ~ 10'. Note the remarkable 
difference between the Ha and UV profiles at (subcritical) galactocentric radii > 5 kpc. 
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Fig. 6. — Median FUV-NUV as a function of galactocentric radius, with individual inner 
(dots) and outer (diamonds) stellar complexes plotted for comparison. We also show the 
FUV-NUV color for model starburst populations of varied age and three different periods 
of continuous star-formation (CSP: 100 Myr, 1 Gyr, and 10 Gyr). Intrinsic and reddened, 
E(B-V) = 0.2, colors are indicated with thick and thin line segments, respectively, although 
the difference is minimal compared to the observed run of FUV-NUV. The < S sas > profile 
(dotted) is duplicated from Fig. 5. 



